Wigner Crystallization in mesoscopic 2D electron systems 
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Wigner crystallization of electrons in a 2D quantum dots is reported. It proceeds in two stages: 
I) via radial ordering of electrons on shells and II) freezing of the inter-shell rotation. The phase 
boundary of the crystal is computed in the whole temperature-density plane, and the influence of 
quantum effects and of the particle number is analyzed. 



In recent years there is growing interest in finite quan- 
tum systems at high density or/ and low temperature. 
In particular, the behavior of a small number of elec- 
trons in quantum dots is actively investigated both ex- 
perimentally HI and theoretically [g, ||. The limiting 
behavior oftwo-dimensional (2D) finite quantum systems 
at zero temperature has been studied by unrestricted 
Hartree-Fock calculations Q which revealed a transition 
from a Fermi liquid to an ordered state called "Wigner 
molecule" . The same crossover at finite temperature has 
been recently demonstrated by fcrmionic path integral 
Monte Carlo jjj.lt has to be expected that further in- 
crease of correlations (increase of the Brueckner param- 
eter r s ) will lead to a still higher ordered quantum state 
resembling the Wigner crystal (WC) 0, 0. 

On the other hand, for finite classical systems, Monte 
Carlo simulations have shown evidence of crystallization 
for sufficiently large values of the coupling parameter T. 
These classical clusters consist of well separated shells 
[B [|, 0), and melting proceeds in two stages: first, ori- 
entational disordering of shells takes place - neighbouring 
shells may rotate relative to each other while retaining 
their internal order. Further growth of thermal fluctua- 
tions leads to shell broadening and overlap - radial melt- 
ing. The temperature of radial melting T r may be many 
orders of magnitude higher than the orientational melt- 
ing temperature T a ||. Large clusters with TV > 100 
have a regular triangular lattice structure and exhibit 
only radial melting. 

Now the question arises, how does the behavior of finite 
electron clusters change at low temperature, i.e. in the 
quantum regime! In this Letter we demonstrate that, 
indeed, Wigner crystallization in 2D quantum electron 
clusters exists and that it is accompanied by two distinct - 
radial and orientiational - ordering transitions too. How- 
ever, in contrast to classical clusters, we observe a new 
melting scenario which is caused by quantum fluctuations 
and exists even at zero temperature ( "cold" melting |u| ) . 
We present a detailed analysis of the two-stage quantum 
melting process and provide numerical data for the phase 
boundaries of both crystal phases, for particle numbers 
in the range N = 10 ... 20. 

Model and characteristic parameters. The theoretical 
analysis of quantum confined electrons at finite temper- 
ature requires the simultaneous account of strong cor- 
relations and quantum effects which excludes e.g. per- 
turbation or mean field methods. We, therefore, use a 



path integral Monte Carlo (PIMC) approach. We con- 
sider a finite unpolarized |ll]] 2D system of N electrons at 
temperature T. The electrons interact via the repulsive 
Coulomb potential and are confined in a harmonic trap of 
strength loq. The system is described by the hamiltonian 
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where m* and are the effective electron mass and 
background dielectric constant, respectively. We use 
the following length and energy scales: Tq, given by 
e /eb^o = m*w 2 rg/2, and E c - the average Coulomb en- 



ergy, 
{r- 



E c = e /efero. After the scaling transformations 
r/ro, E — > E / E c } the hamiltonian takes the form 
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where n is the dimensionless density, n — \f2l\jr\ = 
(ag/ro) 1 / 2 w r^ 1 ^ 2 , a* B is the effective Bohr radius, and 
Iq = h/(m*u>o), is the extension of the ground state 
wave function of noninteracting trapped electrons. (In the 
crystal phase, ro is very close to the mean interparticle 
distance f). The dimensionless temperature is given by 
T = k B T/E c . 

To obtain the configuration and thermodynamic prop- 
erties of clusters of N electrons described by the hamili- 
tonian (||), we performed fermionic PIMC simulations 
using a standard bisection algorithm |D| . The number 
of time-slices M has been varied with n and T accord- 
ing to M = In/T, where I was typically in the range of 
1 ... 10 to achieve an accuracy better than 5% for the 
quantities see below. To obtain the phase bound- 

ary of the Wigner crystal, calculations in a broad range 
of parameter values {n, T, N} were performed. For each 
set of parameters, several independent Monte Carlo runs 
consisting of approximately 10 6 steps were carried out. 

Structure of electron clusters. The simulations yield 
the spatial electron configurations in the trap, examples 
of which are shown in Figs. |l| and |^. One clearly sees 
the formation of shells. Our analysis revealed the same 
shell structures as reported for the corresponding classi- 
cal systems 0. Number of shells and shell occupa- 
tion depend on N, some examples for N — 10 ... 20 are 
given in Tab. [|. Clusters, in which the particle number 
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on the outer shells are multiples of those on the inner 
shells have the highest symmetry which leads to partic- 
ular properties. Examples of these magic numbers are 
N = 10, 12, 16, 19, cf. Tab. |. Let us now discuss the 
influence of quantum effects on the clusters. In contrast 
to classical systems, where the electrons are point-like 
particles, in our case, the wave function of each electron 
has a finite width and may be highly anisotropic, which 
is typical for low temperature, as is most clearly seen in 
Fig. 0a. This peculiar shape results from a superposition 
of N-body correlations, quantum effects and the confine- 
ment potential. Varying the density and temperature, 
the shape changes in a very broad range, which can lead 
to qualitative transitions of the cluster, including cold 
quantum melting, as will be shown below. 

Melting and phase transitions. Wigner crystallization 
is known to occur when the ratio of the Coulomb energy 
to the kinetic energy, E c /Eki n , exceeds a certain thresh- 
old. In the quantum limit, this ratio is just the Brueckner 
parameter r s . On the other hand, in the classical limit 
this ratio goes over to the classical coupling parameter 
T = E c /kBT, In our units (f w r ), 
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Let us first consider the classical case. At sufficiently low 
T (high r), a classical cluster is in the orientationally 
and radially ordered crystal phase. Our simulations re- 
vealed that in a finite system, (in the classical part of the 
phase space), both orientational and radial melting are 
determined by characteristic values of T which, however, 
strongly depend on the particle number N. We find that 
orientational melting occurs at temperatures which vary 
over many orders of magnitude, T* as 3 x 1CP 3 . . . 10~ 12 
and are very sensitive to the shell configuration, as ob- 
served in Ref. ||. Higher temperatures (lower V) are 
found for magic clusters, cf. Tab. Q. The magic clus- 
ter with N — 19 has unusual stability against inter-shell 
rotation, as the ratio of particle numbers on the three 
shells is optimal for the formation of a triangular lattice. 
In contrast, nonmagic clusters have much lower orienta- 
tional melting temperatures (see Tab. |), which is partic- 
ularly striking for A" = 20. 

The situation is completely different at radial melt- 
ing. Here, the critical temperatures for inter-shell par- 
ticle jumps (radial disordering) for magic clusters have 
been found to be approximately 2 times smaller than 
thosefor nonmagic clusters (Tab. |). 

We now investigate the phenomenon of quantum ori- 
entational melting ( OM). As an illustration, consider the 
cluster with A" = 12 particles. Fig. [l] shows three snap- 
shots of particle configurations at low temperature where 
the density increases from(a) to (c), and (b) corresponds 
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FIG. 1: Electron configurations in the x-y plane for N = 12, 
a) in the orientationally ordered state (n = 0.01); b) at ori- 
entational melting (n = 0.014); c) in the orientationally dis- 
ordered/radially ordered state (n = 0.02). Lower right figure 
shows the corresponding angular pair distribution functions 
for electrons on the inner shell. T = 1.0 x 10 -4 . 



tron transitions between two closely lying energy levels. 
At still higher density, Fig. |l](c), the electron wave func- 
tion spreads further, splitting into three maxima, and it 
starts to spread in radial direction. Obviously, OM is 
accompanied by intensification of angular distance fluc- 
tuations, cf. Fig. |l|(d). The spread of the electron wave 
functions smears out the energy barrier for inter-shell ro- 
tation which, in turn, reduces the melting temperature. 
In the limit T — > this leads to so-called "cold" orienta- 
tional melting. 

To allocate the melting temperatures and densities ac- 
curately, we examine the relative mean angular distance 
fluctuations of particlcsfrom different shells 
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where <pi and <f>j are the angular positions of particles on 
shells s\ and S2, respectively, and m Sl , m S2 are the total 
number of particles on the shells. We further consider 
the magnitude of the relative distance fluctuations 
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FIG. 2: Relative angular and radial fluctuations, Eqs. (^[ |^), 
in the vicinity of the orientational (OM) and radial melting 
(RM) for N = 12 and N = 19 versus density. Upper three 
figures show snapshots of the magic cluster with N = 19 in the 
three phases (left n = 0.025, middle n = 0.06, right n = 0.14). 
T = 5.0 x 10~ 4 . Shown error bars are typical for all curves. 



occurs at some critical value of T, see above. On the 
other hand, in the strong quantum limit, there is only 
a weak temperature dependence, so the first method is 
advantagous. To highlight the peculiarities of quantum 
melting, in the following, we concentrate on the density 
dependence of the fluctuations . 

Fig. H shows the n-dependence of the fluctuations 
and u r , at a fixed temperature T = 5 • 10~ 4 , for the 
magic clusters N — 12 and N — 19. The first jump of 
the fluctuations at comparatively low cricitical densities 
n* corresponds to quantum OM. At densities n> n*, the 
angle between particles from neighbouring shells can take 
arbitrary values. There is a significant difference in the 
magnitude of the melting densities for the clusters under 
consideration. As in the classical case, the highest sta- 
bility (highest value n* ) is found for the magic clusters, 
with the maximum value found for N = 19. In par- 
ticular, we obtained the densities of "cold" orientational 
melting by extrapolation of our data to zero temperature. 
The corresponding values for the density and Brueckner 
parameter r£°^ are listed in Tab. 0. 

Let us now proceed in Fig. || to higher densities. We 
observe a clear second jump of the fluctuations u$ and 
u r , which is related to total melting of the WC. It is 
instructive to compare the state of the cluster before and 
after this jump, see middle and right shapshots in Fig. ||, 
respectively. Evidently, when the density is increased, 
the peaks of the wave functions broaden in radial (and 
angular) direction until their width becomes comparable 
with the inter-shell spacing. As a result, the probability 
of inter-shell transfer of electrons grows rapidly, causing a 



sudden increase of the radial fluctuations and the onset of 
radial melting. Our simulations revealed that the jumps 
of and u r are clearly visible along the whole WC phase 
boundary.In the limit of zero temperature, we observe 
"cold" radial melting, the corresponding values for the 
density and rs are presented in Tab. 0. 

Phase boundary of the mesoscopic Wigner crystal. The 
obtained values of the melting temperatures and densities 
are now used to draw the phase boundaries of the radially 
and fully ordered states of mesoscopic clusters. The re- 
sults are summarized in Fig. ^ and Tabs. |, [n]for various 
particle numbers. Consider first the line of radial melting 
"RM". At low densities, n < 0.03, crystallization occurs 
at critical values of T = 1/T*, see the data in Tab. Q. Wc 
found that, for magic (non-magic) clusters, T is above 
(below) the 2D bulk value, = 137, e.g. @. This 
shows that non-magic clusters are more stable against 
radial disordering. The same tendency is observed for 
"cold" melting, where the critical values of r s are system- 
atically lower for non-magic clusters. Interestingly, for all 
comparatively small clusters with 10 < N < 20, the crit- 
ical r s , (Tab. |n|), exceeds the known value of an infinite 
2D system « 37 jL2j . This systematic deviation from 
the critical data of a macroscopic system points to the 
existence of a different scenario of radial melting typical 
for finite systems. For example, for clusters with N = 10 
(19), when approaching the phase boundary, we observe 
increasingly frequent transitions between the shell config- 
urations {2, 8} {3, 7} ({1, 6, 12} *=> {1, 7, 11}). The ex- 
planation is the lowering of the potential barrier for inter- 
shell transitions of individual electrons between these two 
energetically rather close configurations. 

A particularly interesting feature which is missing in 
infinte systems is the existence of the second transition 
characterized by freezing of the inter-shell rotation, cf. 
areas bounded by the lines "OM" in Fig. ||. This highly 
ordered state is confined to the region of low density and 
low temperature. Its phase boundary is highly sensitive 
to the angular symmetry of the cluster. The most sym- 
metric clusters are the magic ones which are essentially 
more stable against orientational melting which is seen 
in their relatively low values of T* and ri , cf. Tabs. | 
and [H| The critical values for the non-magic clusters are 
several orders of magnitude larger, confining their ori- 
cntationally ordered state to much lower densities and 
temperatures. Finally, the last two columns of Tab. || 
contain the maximum temperatures at which the orien- 
tationally and radially ordered crystal phases may exist. 

In summary, we have presented a detailed PIMC- 
analysis of Wigner crystallization of finite electron sys- 
tems confined in a 2D harmonic trap including the critical 
data in the whole density-temperature plane. We have 
shown that there exist two phases characterized by ra- 
dial and radial plus angular ordering, respectively, which 
is essentially different from macroscopic systems. 
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FIG. 3: Phase diagram of the mesoscopic 2D Wigner crys- 
tal. "OM" ( "RM" ) denotes the orientational (radial) melting 
curves for N = 10, 11, 12, 19, 20. Insert shows an enlarged pic- 
ture of the low-density region. Dotted straight lines indicate 
the radial melting transition of a macroscopic classical and 
quantum WC. Brueckner parameter follows from the density 
by r s = 1/n 2 . Shown error bars are typical for all curves. 
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The phase boundaries have been found to be very 
sensitive to the electron number and to the shell sym- 
metry. In contrast to classical clusters, we observed 
"cold" orientational and radial melting which is governed 
by the spread of the electron wave functions in angu- 
lar and radial direction. The predictions of our model 
calculations are expected to be relevant in particular 
for electrons in external fields. Furthermore, they lead 
us to expect Wigner crystallization also in small 2D is- 
lands of electrons [holes] in semiconductor heterostruc- 
tures: for example, in GaAs/AlGaAs systems, crystal- 
lization is predicted for carrier densities below approx- 
imately 10 8 cm -2 [(10 9 . . . 10 10 )cm -2 ] and for tempera- 
tures below 1.6K . . . 5.5AT, for confinement potentials of 
3meV . . . lOmeV. 
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TABLE I: Classical Wigner crystal: critical temperature T* 
and coupling parameter F* corresponding to the orientational 
(o) and radial (r) melting for clusters of size TV in the classical 
limit (n — > 0). Second column contains the shell occupancy, 
starting from the inner shell. Non-magic clusters are typed 
italic. 
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TABLE II: Quantum Wigner crystal: critical density and 
Brueckner parameter r s for orientational (o) and radial (r) 
melting at zero temperature. T™S is the highest possible 
melting temperature, T = n * T. The critical parameters 
n ,r ( s o) and f nax for N = 11 and 20 are estimates @. 
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